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Abstract 

We study how the recent ATLAS and CMS Higgs mass bounds affect the renormalization group 

running of the physical parameters in universal extra dimensions. Using the running of the Higgs 

self-coupling constant, we derive bounds on the cutoff scale of the extra-dimensional theory itself. 

O j We show that the running of physical parameters, such as the fermion masses and the CKM mixing 

Q_[ matrix, is significantly restricted by these bounds. In particular, we find that the running of the 

<D . 

r*j . gauge couplings cannot be sufficient to allow gauge unification at the cutoff scale. 



*Electronic address: Mattias.Blcnnow@mpi-hd.mpg.de 
^Electronic address: mclbcus@kth.se 
^Electronic address: tommy@theophys.kth.se 



Electronic address: hc.zhang@mpi-hd.mpg.de 



I. INTRODUCTION 



Recently, the ATLAS and CMS collaborations presented new bounds on the mass of 
the Standard Model (SM) Higgs boson, excluding values outside the range 115.5 GeV - 
131 GeV and 115 GeV - 127 GeV Q, respectively, at 95 % confidence level. While the 
search for the Higgs boson is the primary goal of the LHC, the experimental collaborations 
are also intensively searching for signs of new physics beyond the SM. Among the most 
popular models describing new physics within the reach of the LHC is the universal extra 
dimensions (UED) model 0]. In this model, all of the SM fields are promoted to a higher- 
dimensional spacetime, giving rise to infinite Kaluza-Klein (KK) towers. The lowest-mass 
KK modes are usually assumed to be located at the TeV scale, and in particular, the lightest 



KK particle could be an interesting dark matter candidate 

An important feature of extra-dimensional models is the impact of the large number of 
KK modes on the renormalization group (RG) running of physical parameters. The RG 
running in extra-dimensional models has previously been investigated, e.g., in Refs. j6-13|. 
It has been shown that the RG evolution changes from the typical logarithmic running in 
four- dimensional models to an effective power-law running at high energies. This means that 
sizable running could take place at relatively low energy scales. In particular, the possibility 
of achieving gauge coupling unification at intermediate energy scales has been discussed 

B0- 

In this paper, we use the RG evolution of the Higgs self-coupling constant in order to 
derive bounds on the minimal five-dimensional UED model, using the recent LHC Higgs 
mass bounds. Previously, results from LHC Higgs searches have been used to constrain five- 
and six-dimensional UED models in Refs. [l4L[15j|, giving the bound R' 1 > 700 GeV for the 
minimal five-dimensional UED model. This corresponds to a n up per bound ran < 500 GeV 
on the Higgs boson mass from electroweak precision data [l6j. The RG running of the 



Higgs se 
Refs. 
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f-coupling constant has also recently been used to constrain new physics models in 



19]. 



In addition, we discuss the running of fermion masses and mixing parameters in the 
UED model, taking the new bounds into account. These fundamental physical parameters 
are crucial for building new physics models, as well as testing the feasibility of theories 
beyond the SM. In fact, since the values of these parameters are not predicted by the SM, 
new physics, which is usually located at some very high energy scale, is needed in order to 
gain insight into their origin. Thus, we provide values for the fermion masses and mixing 
parameters at the cutoff scale of the UED model. 

The rest of this work is organized as follows: In Sec. [Til we discuss general features of 
renormalization group running in extra-dimensional theories. Next, in Sec. we discuss 
the running of the Higgs self-coupling constant and the resulting bounds on the UED model 
from the LHC Higgs mass bounds. Then, in Sec. IIVI we show the running of the gauge 
coupling constants and demonstrate that gauge unification cannot be achieved within the 



2 



UED model. In Sees. |V] and IVIl we give the RG evolution of the fermion masses and the 
CKM matrix parameters, respectively. Finally, in Sec. IVII[ we summarize and discuss our 
results. In addition, in Appendix |AJ we provide the one-loop beta functions that are relevant 
for our work. 



II. RENORMALIZATION GROUP RUNNING IN EXTRA DIMENSIONS 

A general feature of quantum field theories with extra spatial dimensions is that they are 
non-renormalizable. However, as pointed out in Ref. 7|> such models could preserve renor- 
malizability if they are truncated at a certain energy scale (i.e., the number of KK modes is 
finite). In such a situation, physical quantities are subject to a power-law running behavior, 
in contrast to the typical logarithmic running in ordinary four- dimensional theories. This 
power-law running may significantly change the running physical parameters, such as the 
gauge couplings, the quark and lepton masses and mixing angles, and the Higgs mass. 

In general, the beta function of a parameter P in an extra-dimensional model can be 
expressed as 

HP 

167r 2 -— = (3 P + s(3 P , (1) 
dm/i 

where (3p denotes the SM beta function, while 0p corresponds to the contributions from the 
KK modes at any single KK level to the total beta function. Here, it is assumed that the 
particle content at each non-zero KK level is the same, except for the particle masses. This 
is the case in many models, and in particular in the UED model. The scale parameter s 
is defined as s = [p/fJ>o\, where \_x\ is the largest integer smaller than x, and /xo = -R -1 is 
the inverse radius of the extra dimensions, i.e., s counts the number of KK levels below the 
energy scale /i. At energy scales below /io, i.e., below the mass of the lowest KK excitations, 
the (3p term can be ignored, whereas for /i ^> /i , many KK modes are excited and their 
contributions change the scale-dependence of the physical parameters from logarithmic to 
power-law. The relevant one-loop beta functions for the five- dimensional UED model can 
be found in Appendix [A] 

As mentioned above, in order to make the theory renormalizable, an explicit cutoff scale 
A has to be introduced. From this point of view, the UED model is an effective description 
at low-energy scales, which is replaced by a renormalizable theory above the cutoff scale. 
In the UED model, A is usually taken to be the energy scale where the gauge couplings 



become non-perturbative 20j, but it could also be related to a unification scale for the 



gauge couplings jsl, 0). In this work, we will apply the LHC bounds on the Higgs mass to 
the running behavior of the Higgs self-coupling constant in order to put bounds on A. 

In the following sections, we perform a numerical analysis of the running physical pa- 
rameters. In our computations, we make use of the full one-loop RGEs without any fur- 
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ther approximations. The input values for the physical parameters, at the energy scale 
M z = 91.2 GeV, are taken from Ref. [21]. 



III. RUNNING OF THE HIGGS SELF-COUPLING CONSTANT 

The running of the Higgs self-coupling constant A in the UED model is given in Eq. ( 1A1I) . 
Now, the RG evolution of A can be used to constrain the UED model. In particular, for 
different initial conditions, A may approach the triviality limit (A diverges), or the vacuum 
stability limit (A becomes negative, i.e., the Higgs potential becomes unstable). Note that 
a negative A does not necessarily mean that the model is invalid, since the electroweak 
vacuum might be metastable. In other words, we could live in an unstable vacuum, while 
the lifetime of the electroweak vacuum is longer than the age of the Universe. However, 
the metastability limit relies on the fastest process conceived for the transition to the true 
vacuum, and any faster process occurring once anywhere in the Universe would reduce or 



eliminate the metastability region. For a general discussion on this topic, see, e.g., Ref. 22]. 
Using the relation A = m 2 H /v 2 , this allows us to constrain the parameter space of the UED 
model using the LHC Higgs mass bounds, by requiring that neither the triviality nor the 
vacuum stability limit is reached below the cutoff scale. In the low Higgs mass region 
allowed by the new LHC bounds, only the vacuum stability condition is important for the 
UED model. 

In Fig. [U we present the bounds on the cutoff scale A from the requirement of vacuum 
stability. We show upper bounds on the product AR, which counts the number of KK levels 
below the cutoff scale, as a function of Rr 1 . The results depend on the value of the top 
quark mass m t , which is only known to an accuracy of a few GeV, and therefore, we present 
our results for m t in the range 170.9 GeV - 173.3 GeV. The weakest bounds are obtained in 
the phenomenologically interesting range around 1 TeV. We observe that the global upper 
limit on the number of KK modes in the model is five only, constraining the validity range 
of the extra-dimensional description significantly. As we have mentioned, the metastability 
of the SM vacuum could result in a smaller lower bound on the Higgs mass, which is the 
situation in the UED model. In a semi-classical approximation for A > —0.05, the SM Higgs 
vacuum may still be stable against the age of the Universe {23} , and thus not in conflict with 
experimental observations. This will actually increase slightly the upper bound on AR in 



Fig. [TJ See Ref. [19| for a detailed numerical analysis on this issue. It should also be noted 
that our results rely on the one-loop beta functions, and would be slightly changed by taking 
higher-order contributions into account. Nevertheless, we expect our main conclusions to 
remain valid at higher order. 
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FIG. 1: Upper bounds on AR as a function of -R -1 from the vacuum stability condition for the 
Higgs self-coupling constant. The bands show the variation of the bound with the top mass in the 
range 170.9 GeV - 173.3 GeV, where the strongest bounds are obtained for the largest value for 
m t . 



IV. RUNNING OF THE GAUGE COUPLINGS 

Next, the running of the gauge couplings in the UED model is given in Eq. (1A4|) . Solving 
this equation, we obtain 



1 



1 



9t(M) 9f(M z ) 8tt 
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s In 



Ho 



— In s\ 



(2) 



between the n-th and (n + l)-th thresholds. The second term in Eq. (j2J) corresponds to 
the SM contributions and the last term to the corrections from the KK modes. Note that 
the expression in the last parenthesis is always positive, e.g., for \ij \i§ = 5,10,40 one has 
sln(/i//! ) — Ins! ~ 3,8,37. In the limit of large s, i.e., ji ^> /i , the evolution of is 
dominated by the contributions from KK excitations. Since the coefficients hi and hi are in 
general not the same, the impact on the RG running from the KK modes is different from 
that from the SM particles. In particular, the sign of 62 is opposite to that of 62, indicating 
that Q2 tends to increase at higher energy scales. In addition, h\ = 27/2 is larger than the 
other two coefficients, which leads to a fast running behavior of g\ at high energies. 

By solving the RGEs for the gauge couplings, we obtain the running of the gi, which is 



5 



shown in Fig. [2] for i?" 1 = 1 TeV. The most interesting feature of this result is that the 
region where the coupling constants would approximately unify is ruled out by the vacuum 
stability criterion. In fact, it turns out that this is a general statement and there is a no-go 
scenario for gauge unification below the cutoff scale A in the UED model. From Eq. (j2J), it 
follows that 



DyQi) = Dy(M z ) - ^ j&y In + b {i \s In \T-J - hislj I EE £>y(M Z ) - A^/l), 

„ (3) 

where = 1/gf — l/g?, fey = bi — bj, and &y = h — bj. By comparing Ay (A) with Dij(Mz), 
we can observe from Fig. [3] that the ratio never reaches (or is even close to) one. In this 
figure, we have taken the upper limit of the cutoff scale from Fig. [1] for m# = 130 GeV, 
which is equivalent to using the global upper limit on A. Thus, Dij(/i) can never become zero 
below the cutoff scale, meaning that the gauge couplings gi and gj will not unify while the 
extra-dimensional theory is valid. Note that, while higher-order corrections or allowing for 
a broader uncertainty range in the input parameters could change the actual values slightly, 
the ratio is quite far away from one and this conclusion should therefore be robust to such 
details. 




\x [GeV] 

FIG. 2: The RG evolution of the three gauge couplings as functions of the energy scale in the UED 
model, with R^ 1 = 1 TeV. The gray-shaded area is ruled out by the vacuum stability criterion. 
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FIG. 3: The ratio Aij(A)/Dij(Mz), where A is given by the global upper limit from Fig. [TJ The 
solid, dashed, and dotted curves correspond to {i,j} = {1,2}, {1,3}, and {2,3}, respectively. A 
ratio of one or higher would indicate that the gauge couplings g^ and gj unify below A, which this 
figure shows is impossible given the current bounds on the Higgs mass. 

V. RUNNING OF THE FERMION MASSES 

As mentioned before, the evolution of the fermion masses is important for constructing 
new physics models. Therefore, in this section, we give a detailed discussion on the running 
fermion masses at the cutoff scale for various radii R. 

The running quark masses can be obtained from Eq. (IA5p as 

rhf = [Re(F u )f f + a u + sa u ] m f , (4) 
where rhf = 167r 2 dm//dln/i, / = u,c,t, and F u is given by 

F u = ^(D u -VD d V^)(l + s) . (5) 
For the down-type quarks, a similar relation holds, with 

F d = l(V^D u V -D d ){l + s) . (6) 
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Here, V denotes the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and 



D u = diag(yl,y 2 c ,y 2 ) 
D d = diag{y 2 d ,y 2 ,y 2 ) 



(7) 
(8) 



where the y 2 are the eigenvalues of the matrix YjYf. We adopt the standard parametrization 
of the CKM matrix, in which V is parametrized by three mixing angles and one CP-violating 
phase, viz., 



V 



C12C13 



S12C13 



sue 



-S12C23 - ci 2 s 2 3Si 3 e l5 C12C23 - si2S 2 3Si 3 e 1<5 S23C13 

\^ S12S23 — Ci2C23Si3e l5 — C12S23 — Si2C23Si 3 e l5 C23C13 J 



(9) 



with Cij = cos9ij and Sij = sin% = {1,2}, {1,3}, {2,3}). Note that, above the 

electroweak symmetry breaking scale, the unbroken gauge symmetry forbids quark and 
lepton masses. The actual meaning of a fermion mass rrif in this energy region is a measure 
of the corresponding non-trivial Yukawa coupling eigenvalue yj. We adopt the definition 
rrif = yfV, above the electroweak scale, where v ~ 174 GeV is the vacuum expectation value 
of the Higgs field in the SM. 

In view of the hierarchical spectrum of the quark masses, i.e., m t ^> m^ c ^> m u ^ s , one 
can neglect all of the Yukawa couplings except for y t , and in this approximation, the running 
of the quark masses is given by the equations 



rn u ~ (a u + sa u ) m u , 
m c ~ (a u + sa u ) m c , 

^ (Vt + s Vt ) + a u + sa u 



m t 



(10) 
(11) 
(12) 



and 



m d 



-( 

2 v 

-( 



2 2 2 
S 13 C 12 C 23 



2 2 2 
S 12 S 13 C 23 



s 2 12 s 2 23 - 2s 12 s 13 S23C 12 C23Cs) (y 2 + sy 2 ) +a d + sa d 
s 23 c i2 + 2S12S13S23C12C23Q) (y 2 + sy 2 ) +a d + sa d 



m d 



2 C 23 c is {Vt + s Vt) + «rf + sa d 



m b . 



(13) 
(14) 
(15) 



The RG evolution equations for m u and m c are similar to each other and are both governed 
by the flavor-diagonal part a/, whereas for the top quark mass, contributions from y t should 
be taken into account, which slightly changes the running of m t . In the down-type quark 
sector, the two light quarks d and s also receive similar RG corrections, since the flavor 
non-trivial parts [i.e., the yt terms in Eqs. ffl3|) and (fT4"j) ] are suppressed by the CKM mixing 
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angles. Since the «/ parameters are negative, we expect the quark masses to decrease with 
the energy scale. 

As for the charged leptons, we can safely ignore the Yukawa corrections due to the 
smallness of their masses, and we obtain 

rhi = (at + sa^rrii , (16) 

for % = e,/j,,T, where a? is the flavor-diagonal part of the right-hand side of Eq. flA9[) . 
Therefore, the charged-lepton masses are essentially rescaled by a common factor at high- 
energy scales. However, this factor is larger than in the SM, due to the scale parameter s. 
Furthermore, in contrast to the quark sector, the flavor- diagonal part is positive, due to 
the lack of #3 corrections (i.e., leptons do not participate in the strong interactions), which 
leads to larger values for the charged-lepton masses at higher energies. 

In order to numerically show the RG evolution of the fermion masses, we define the ratios 
Rf = rrif(n) /rrif(Mz) reflecting the RG corrections at the scale /x. The scale-dependence of 
m t , 777,6, an d m T is illustrated in Fig. HJ for R" 1 = 1 TeV and m# = 125 GeV. As expected, 
the running quark masses decrease with the energy scale. We observe that the quark masses 
at the cutoff scale are reduced by about 15 % - 20 % (c.f., Table fl] for detailed numbers). 
Therefore, it seems impossible to achieve a reasonably good unification of Yukawa couplings, 
due to the bounds on A. As for the lepton sector, the charged-lepton masses run to larger 
values at the cutoff scale A, e.g., a 10 % increase of the tau mass can be observed from 
the plot. It is also interesting to point out that the values for the running charged-lepton 



masses 



24j are maximal at /x = 30 TeV - 40 TeV (or equivalently n = 30 - 40), a region 
which is not allowed by the new Higgs mass bounds. It should be stressed that the running 
of the charged-fermion masses is not very sensitive to the specific value for the Higgs mass, 
since A does not enter the beta functions for Yf at one-loop level. As a reference for model 
building, we list in Table U the running quark and charged-lepton masses at the cutoff scale 
A for R~ l = 500 GeV, 1 TeV, and 10 TeV. These numerical values are consistent with Fig. HJ 
and in general, the evolution of the masses is relatively small due to the strong constraints 
on A. We hope that the table could be useful for building possible extra-dimensional models 
within the UED framework. 

Finally, we investigate the running neutrino masses with the new Higgs mass bounds in 
the UED model. By analytically diagonalizing the beta function for k, which is given in 
Eq. flA6j) . we arrive at very compact expressions for the evolution of the neutrino masses 

m - (ot K + sa R )mi , (17) 

where a K is the diagonal part of the right-hand side of Eq. (IA10I) and we have omitted the 
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FIG. 4: The RG evolution of the mass ratios Rf for the top quark, the bottom quark, and the tau 
lepton in the UED model, for R~ 1 = 1 TeV and mu = 125 GeV. The gray-shaded area is ruled 
out by the vacuum stability criterion. 

charged-lepton Yukawa couplings. An approximate solution to this equation is given by 



Therefore, the RG running of the neutrino masses is only sensitive to a K , independently of 
the neutrino mass spectrum and the mixing parameters. Similarly to the charged fermions, 
we define the ratios Ri = rrii(fi) / rrii(Mz) for the neutrino masses, and illustrate the evolution 
of Ri in Fig. [5] for m# in the range 115 GeV - 130 GeV. An important feature of the running, 
which can be seen from the plot, is that, due to the stability bounds, Ri cannot reach large 
values below A. Furthermore, the running of the neutrino masses does indeed depend on the 
Higgs mass tuh, since the effective neutrino coupling matrix k receives one- loop corrections 
from the quartic Higgs interaction. 

VI. RUNNING OF THE CKM MIXING MATRIX 

The CKM mixing matrix stems from the mismatch between the diagonalization of the 
Yukawa matrices Y u and Yd, and the running of the CKM matrix is not sensitive to the 
flavor- diagonal parts in the beta functions for Yf. Explicitly, one could insert the CKM 




(18) 
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177 , [Me VI 


2.82 


2.4 


2.3 


2.1 


m u [MeV] 


1.38 


1.2 


1.1 


1.0 


m T [GeV] 


1.746 


1.9 


1.9 


1.9 


[MeV] 


102.7 


110 


110 


110 


m e [MeV] 


0.4866 


0.53 


0.53 


0.52 



TABLE I: Fermion masses at the cutoff scale A for R~ l = 500 GeV, i?" 1 = 1 TeV, and R' 1 = 
10 TeV, respectively. The input values for the fermion masses at the energy scale [i = M z are 
listed in the left column for reference. 

matrix into Eq. (IA5[) . and obtain the individual beta functions for the CKM mixing angles 

3 

^12 = ~2 (Vt + S Vt ) c 12 [(s? 3 C2 3 - S 2 23 ) Sia + 2S23S13C12C23C5] , (19) 
023 = \ (Vt + SVt) S23C23 , (20) 

3 

^13 = g + S2/ *) Sl 3 C 13 C 23 > ( 21 ) 

as well as 5 ~ at leading order. 

Note that current experiments indicate that all the quark mixing angles are relatively 
small and the CKM matrix takes a nearly diagonal form. Thus, the beta function for 9 12 
[the right-hand side of Eq. (!l~9~]) ] is strongly suppressed by sines of the mixing angles, implying 
that 9 12 is stable against radiative corrections. The mixing angles 6*23 and #13 may receive 
visible RG corrections, and they increase with the energy scale. In fact, using Eqs. (120]) and 
(121]) . it holds that #23 and 6*13 are related to each other by sm26>i 3 = C tan #23, where C is a 
constant. 

In analogy with the mass ratios, we define the ratios Aij(fi) = 9ij([i)/9ij(Mz) char- 
acterizing the running behavior of the quark mixing angles. The evolution of the Aij is 
demonstrated in Fig. [6] for R^ 1 = 1 TeV. In agreement with our analytical results, 9 12 is 
rather stable, whereas #23 and 9 13 could increase by about 5 %. 
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FIG. 5: The RG evolution of the neutrino mass ratios B4 in the UED model, for ran in the range 
115 GeV - 130 GeV. The gray-shaded area is ruled out by the vacuum stability criterion. 

VII. SUMMARY AND CONCLUSIONS 

In this work, we have studied the RG running of physical parameters in the five- 
dimensional UED model. In particular, we have investigated the impact of the recent 
ATLAS and CMS Higgs mass bounds on the cutoff scale of the extra-dimensional model. 
These bounds come from the criterion of not reaching the Higgs vacuum instability limit 
below the cutoff scale. We have found that the five-dimensional UED model can be valid 
only at most up to the fifth KK level, significantly constraining the higher-dimensional 
description. 

Using this new result, we have shown that it is generally not possible to achieve gauge 
coupling unification at the cutoff scale in the UED model. Furthermore, we have studied the 
RG running of the quark and lepton masses and mixing parameters in the UED model. We 
have found that, while the running at high-energy scales shows interesting features, these 
regions are excluded by the new bounds on the model. In particular, the regions of large 
power-law running is excluded. 

Our results demonstrate that the LHC searches for the SM Higgs boson can have impor- 
tant consequences also for models of physics beyond the SM. As the bounds on the Higgs 
mass become even stronger, the global limit on the cutoff scale could be decreased suffi- 
ciently to allow for only three KK levels in the model. We emphasize that this is only an 
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FIG. 6: The RG evolution of the ratios A^ in the UED model with R = 1 TeV. The curves of 
A23 and A13 overlap with each other. The gray-shaded area is ruled out by the vacuum stability 
criterion. 

upper limit on the cutoff scale, and that new physics that changes the evolution of the Higgs 
self-coupling constant sufficiently to avoid the vacuum instability has to be introduced below 
this scale. 

Although we have considered only the five-dimensional UED model, we expect that even 
stronger constraints can be derived for six-dimensional models, since a higher density of 
states gives rise to an even faster running in such models. 
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Appendix A: One-loop beta functions in the UED model 

The running of the Higgs self-coupling constant is given by 



16tt 2 



dA 



P\ + s/3a • 



(Al) 



dln/i 



13 



Here, we have defined the scale parameter s = \_fi/fio\ , where [^J is the closest integer below 
x. The SM contribution reads 25 



fix = 6A 2 - A (3gl + 9gf) + Qrf + Zg\gl + ^) 

3 {YXT + 3 <YjY d y + (Y}Y t 



+ 4AT - 8tr 



(A2) 



where T 
are 



10| 



tr \ 3Y^Y U + SYJYj + Y^Y^j . In addition, the extra-dimensional contributions 



fix = 6A 2 - A (3gj + 9g 2 2 ) + (2g{ + 4g 2 g 2 2 + Qg*) 



8AT - 16tr 



3W 2 + 3(*) 2 + (y/y / 



(A3) 



Next, the RGEs for the gauge couplings are given by 

2 dg { 



16tt' 



dln/i 



6j + sk g 



3 



(A4) 



where (61,62,63) = (41/6,-19/6,-7) and (61,62,63) = (27/2,7/6,-5/2) (l0|. 

Finally, the one-loop RGEs for the Yukawa coupling matrices Yf (/ = tt, d, £) and the 
neutrino mass operator k can be expressed in a general form as 



16vr 
16tt 2 



2 dY, 



dln/i 
dln/i 



Pf + Pf = Pf + sa f Y f + sY^, , 

/3 K + /3 K = fi K + soL^K + sKiV re + siVjft , 



(A5) 
(A6) 



where the SM beta functions are 



26- 



fiv 

fid 
Pi 
p. 



17 
12' 



Y d {-\Y^Y U + ^Y^ - ±rf - ^ 2 - 8g 2 3 + t) 

^G r ^-f^-^ 2+T ) ' 

k + (A - 3^ + 2T) « . 



k y;y> 



(A7) 
(A8) 
(A9) 
(A10) 
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The contributions from the KK excitations are given by [ll|, 



and 



N u = lY^Y u -^YjY d , (All) 

N d = -\YiY u + \YlY di (A12) 

N £ = \y}Y t , (A13) 

N K = ~\y}Yu (A14) 



101 2 15 2 28 2 

= —^91 - - y# 3 + 2T > ( A15 ) 

& d = -^gt-fgi-fg! + 2T, (A16) 

= -f^ 2 -^ 2 2 + 2T, (A17) 

a K = - l -gl - l -^g 2 2 + 4T + A . (A18) 



[1] ATLAS, ATLAS-CONF-1010.55222011-163 (2011). 
[2] CMS, CMS-PAS-HIG- 11-032 (2011). 

[3] T. Appelquist, H.-C. Cheng, and B. A. Dobrescu, Phys. Rev. D64, 035002 (2001), hep- 
ph/0012100. 

[4] G. Servant and T. M. P. Tait, Nucl. Phys. B650, 391 (2003), hep-ph/0206071. 
[5] H.-C. Cheng, J. L. Feng, and K. T. Matchev, Phys. Rev. Lett. 89, 211301 (2002), hep- 
ph/0207125. 

[6] K. R. Dienes, E. Dudas, and T. Gherghetta, Phys. Lett. B436, 55 (1998), hep-ph/9803466. 
[7] K. R. Dienes, E. Dudas, and T. Gherghetta, Nucl. Phys. B537, 47 (1999), hep-ph/9806292. 
[8] S. A. Abel and S. F. King, Phys. Rev. D59, 095010 (1999), hep-ph/9809467. 
[9] G. Bhattacharyya, S. Goswami, and A. Raychaudhuri, Phys. Rev. D66, 033008 (2002), hep- 
ph/0202147. 

[10] G. Bhattacharyya, A. Datta, S. K. Majee, and A. Raychaudhuri, Nucl. Phys. B760, 117 

(2007), hep-ph/0608208. 
[11] A. S. Cornell and L.-X. Liu, Phys. Rev. D83, 033005 (2011), 1010.5522. 
[12] A. S. Cornell and L.-X. Liu, Phys. Rev. D84, 036002 (2011), 1105.1132. 
[13] M. Blennow, H. Melbeus, T. Ohlsson, and H. Zhang, JHEP 04, 052 (2011), 1101.2585. 
[14] K. Nishiwaki, K.-y. Oda, N. Okuda, and R. Watanabe, Phys. Lett. B707, 506 (2012), 

1108.1764. 



15 



[15] T. Kakuda, K. Nishiwaki, K.-y. Oda, N. Okuda, and R. Watanabe (2012), 1202.6231. 

[16] M. Baak, M. Goebel, J. Haller, A. Hoecker, D. Ludwig, et al. (2011), 1107.0975. 

[17] M. Holthausen, K. S. Lim, and M. Lindner, JHEP 02, 037 (2012), 1112.2415. 

[18] I. Masina and A. Notari (2011), 1112.2659. 

[19] J. Elias-Miro, J. R. Espinosa, G. F. Giudice, G. Isidori, A. Riotto, et al., Phys. Lett. B709, 

222 (2012), 1112.3022. 

[20] D. Hooper and S. Profumo, Phys. Rept. 453, 29 (2007), hep-ph/0701197. 

[21] Z.-z. Xing, H. Zhang, and S. Zhou (2011), 1112.3112. 

[22] M. Sher, Phys. Rept. 179, 273 (1989). 

[23] G. Isidori, G. Ridolfi, and A. Strumia, Nucl. Phys. B609, 387 (2001), hep-ph/0104016. 

[24] Z.-z. Xing, H. Zhang, and S. Zhou, Phys. Rev. D77, 113016 (2008), 0712.1419. 

[25] T. P. Cheng, E. Eichten, and L.-F. Li, Phys. Rev. D9, 2259 (1974). 

[26] K. S. Babu, C. N. Leung, and J. T. Pantaleone, Phys. Lett. B319, 191 (1993), hep-ph/9309223. 

[27] P. H. Chankowski and Z. Pluciennik, Phys. Lett. B316, 312 (1993), hep-ph/9306333. 

[28] S. Antusch, M. Drees, J. Kersten, M. Lindner, and M. Ratz, Phys. Lett. B525, 130 (2002), 
hep-ph/0110366. 



16 



